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In this work we study symmetry fractionalization of vison excitations in topological Z 2 spin 
liquids. We show that in the presence of the full SO(3) spin rotational symmetry and if there is an 
odd nnmber of spin-| per unit cell, the symmetry fractionalization of visons is completely fixed. On 
the other hand, visons can have different classes of symmetry fractionalization if the spin rotational 
symmetry is reduced. As a concrete example, we show that visons in the Balents-Fisher-Girvin Z 2 
spin liquid have crystal symmetry fractionalization classes which are not allowed in SO (3) symmetric 
spin liquids, due to the reduced spin rotational symmetry. 


Global symmetries, including the spin rotational sym¬ 
metry, the time-reversal symmetry and the crystal sym¬ 
metries, play important roles in the study of topological 
spin liquids m- Not only are they useful in no-go the¬ 
orems that guarantee the existence of topological orders 
and constrain their properties when certain symmetries 
are unbroken IMD], but they also further enrich m the 
intrinsic topological orders of the spin liquids [MS], which 
are characterized by anyon excitations with fractional 
braiding statistics. Particularly, in the resulting symme¬ 
try enriched topological (SET) order the anyon excita¬ 
tions can exhibit symmetry fractionalization, where the 
anyons carry fractional quantum numbers of the symme¬ 
try group Examples of symmetry fractionaliza¬ 

tion include fractionally charged quasiholes in fractional 
quantum Hall states US], and spinons carrying a half¬ 
integer spin in spin liquids (2 HD] . 

Symmetry fractionalization is crucial in classifying, de¬ 
tecting and modeling the topological order of Z 2 spin 
liquids. First, it refines the classification of topologi¬ 
cal orders in spin liquids: although the anyon content 
is identical, different patterns of symmetry fractionaliza¬ 
tion belong to different quantum phases [I1HI1H2- Sec¬ 
ond, fractional quantum numbers carried by anyons can 
provide smoking-gun experimental signatures to detect 
the anyons and the topological orders |23H25]- Lastly, 
symmetry fractionalization plays a key role in the con¬ 
struction of mean field miiiD] and variational wave func¬ 
tions HUH?] for symmetric spin liquids. 

In a Z 2 spin liquid, symmetry fractionalization is fully 
characterlized by the symmetry quantum numbers car¬ 
ried by spinons (either the bosonic or fermionic ones) 
and visons, respectively |12j . The symmetry fraction¬ 
alization of spinons has been extensively studied using 
parton constructions [2T1 HD] [28D32] and tensor network 
states m- However in all these studies it turns out 
that the symmetry fractionalization of visons is com¬ 


pletely determined by the background spinons on the 
lattice. Recently it has been shown that certain SET 
phases where spinons and visons both have nontrivial 
symmetry fractionalizations are anomalous |33H36] . im¬ 
plying that these phases cannot possibly exist in any two- 
dimensional (2D) systems and in fact can only be realized 
on the surface of three-dimensional (3D) symmetry pro¬ 
tected topological (SPT) phases [37], such as 3D topolog¬ 
ical insulator/superconductors (TSC) [35] and topologi¬ 
cal cystalline insulators (TCI) |3D|- However, although 
there is a general framework to detect anomalous SET 
orders for onsite unitary symmetries [131140] , for SET or¬ 
ders with time-reversal symmety and crystal symmetries 
a general method of detecting anomaly is still lacking. 
Since evidences for gapped Z 2 spin liquid phases with 
full spin rotational, time reversal and crystal symmetries 
have been found in numerical studies of spin-^ Heisen¬ 
berg models on the kagome lattices [411143] (see however 
Ref. I44II , and also in a recent experiment in the spin- ^ 
kagome antiferromagnet herbertsmithite [45] , it is an im¬ 
portant theoretical question to study the classification of 
(anomaly-free) symmetry fractionalization in fully sym¬ 
metric gapped Z 2 spin liquids. 

In this work we study symmetry fractionalization for 
vison excitations in a Z 2 spin liquid with spin rotational, 
time-reversal and crystal symmetries. Generalizing an 
elegant method to detect anomalies recently proposed 
by Hermele and Chen [36] . we show that on kagome or 
square lattices, with full SO (3) spin rotational symmetry 
and under the assumption that the spinon carries a half¬ 
integer spin, all symmetry fractionalization classes for 
visons are uniquely fixed, except the commutation rela¬ 
tion between the two unit translations. In fact, any SET 
orders where visons exhibit other nontrivial symmetry 
fractionalization are anomalous. Our conclusion crucially 
depends on the full SO (3) spin rotational symmetry, and 
we demonstrate using the Balents-Fisher-Girvin(BFG) 
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model |46| that the visons can exhibit more complicated 
fractionalization patterns once the spin rotational sym¬ 
metry is reduced to 0(2). 

Symmetry fractionalization of vison in S0(3) symmet¬ 
ric spin liquids. We consider Z 2 spin liquids enriched 
by the symmetry group SO (3) x x Ggpace, denoting 
the spin rotational symmetry, the time-reversal symme¬ 
try and space group of the lattice, respectively. We will 
focus on the kagome lattice and treat the more compli¬ 
cated case of square lattice in Sec.|IV]of the Supplemental 
Material. Notice that for systems consist of half-integer 
spins, the symmetry group should still be considered 
S0(3) instead of SU(2) because local excitations, such 
as the magnons, always carry integer spins. We will also 
assume that there are even number of sites in the system, 
so the ground state can be a spin singlet. 

Before we discuss the fractionalization of space group 
and time-reversal symmetries, we first determine the 
S0(3) spin quantum numbers of the anyons. A Z 2 spin 
liquid has four types of topologically distinct quasipar¬ 
ticle excitations: the trivial excitation 1, the bosonic 
spinon e and the vison m which are mutually semionic, 
and the fermionic spinon e = exm. Here we assume that 
the e anyon in the Z 2 spin liquid carries a half-integer spin 
projective representation of the SO(3) symmetry group. 
With this assumption one can further show that the time- 
reversal symmetry forces the vison m to have an integer 
spin, otherwise the system must have a nontrivial Hall re¬ 
sponse associated with the charge, which necessarily 
breaks the time-reversal symmetry [33l . On a lattice 

with an odd number of spin-f per unit cell (which in¬ 
cludes most of the gapped spin liquids found in numerical 
simulations so far [HHISI EHHil], and the Lieb-Schultz- 
Mattis-Oshikawa-Hastings theorem PEiinni guarantees 
that any gapped symmetric states found in this sort of 
system must be topologically ordered), it can be shown 
that [HI Hn]: (a) our assumption of e carrying a half¬ 
integer spin is automatically guaranteed, and (b) the vi¬ 
son sees a tt-Aux when moving around a unit cell, i. e. the 
vison transforms projectively under translation symme¬ 
tries: T 1 T 2 = —T 2 T 1 where T 1.2 are translations by the 
two basis vectors of the Bravais lattice, because when 
vison moves around a unit cell it braids with an odd 
number of spinons inside the unit cell, and each braiding 
gives a —1 Berry phase. 

Our argument of fixing vison’s symmetry fractional¬ 
ization is based on the flux-fusion anomaly test recently 
proposed by Hermele and Chen [3^. We will consider 
systems on a disk or an infinite plane. In this test, 
we adiabatically insert fluxes of the U(l) global sym¬ 
metry of spin rotations with respect to S'^. A flux 
can be inserted by applying the unitary transformation 
U = where r labels lattice sites and 9^ de¬ 

notes the polar angle of r with respect to a branch cut, 
as shown in Fig. Because e and e anyons carry half¬ 
integer spins, a 27r flux carries the topological charge of 



FIG. 1. Symmetry actions on fluxes. The cross labels the 
position of fluxes, and the wiggly lines mark the positions of 
the branch-cuts, (a) An antiunitary symmetry X maps a flux 
to itself, (b) A unitary symmetry X maps two fluxes with 
opposite f to each other. 

a vison, ^12tt ~ rn |131136] . If the symmetry fractionaliza¬ 
tion of m can be deduced from a physical invariant of 27r 
fluxes, which can be extended to a well-defined quantity 
for the whole adiabatic flux-insertion process, based on 
continuity we can immediately relate the symmetry frac¬ 
tionalization of ^ 2 tt ~ m to that of the vacuum sector 1, 
which is always trivial. |53| 

We will now expand on this general method, focusing 
on an order 2 symmetry operation X. As we will see 
later this is sufficient for our purpose. 

For an anti-unitary Z 2 symmetry operation A, we im¬ 
plement X such that a flux at position r is mapped 
to itself under A: 

A : VI^{y) H0(r), (1) 

as shown in Fig. |l(a)[ and thus one can define the lo¬ 
cal action of A. Two topologically distinct ways A can 
act are given by A^ = ±1 on the flux, corresponding to 
the flux being a Kramers singlet/doublet under A. For 
0 < (/ < 27r, mathematically this is the well-known fact 
that the second group cohomology with U(l) coefficient 
i/^(Zf-, U(l)) = Z 2 ( Zf- denotes the group generated by 
A), which classifies projective representations of Z^ [37] . 
In the end of flux insertion when a vison is created, due 
to the Z 2 fusion rule m x m = 1 the local A action is clas¬ 
sified by i7^(Zf-,Z2) = Z 2 [T31 [T3] . distinguished again 
by the local A^ value. Therefore the two classifications 
coincide, and the local A^ value is well defined for the 
whole flux-insertion process. We therefore conclude that 
m must have A^ = +1. 

For an unitary Z 2 space group symmetry operation 
A, we can not use the same argument when the flux 
is mapped to itself under A. The reason is that for 
Q < (j) < 27r, the local value of A^ is no longer a 
physical invariant because one can redefine the local ac¬ 
tion of A by an arbitrary U(l) phase, consistent with 
i7^(Z^,U(l)) = TL\. But for a vison, A^ = ±1 are still 
topologically distinct due to the constraint from the Z 2 
fusion rule in accordance with i7^(Z^,Z2) = Z 2 . There¬ 
fore for an unitary symmetry operation we need to con¬ 
sider a different setup where A maps a flux located at 
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FIG. 2. (a) Symmetry operations on a kagome lattice. Ti 
and T 2 denote translations by the two basis vectors of the 
Bravais lattice; and a denote mirror reflections with re¬ 
spect to the axes shown by the dashed line; Ip and Is de¬ 
note plaquette-centered and site-centered inversions (i.e. 180- 
degree rotation), respectively, (b) Duality mapping from the 
BFG model [46] to the quantnm dimer model. The spins on 
the kagome lattice sites are mapped to dimers on the bonds 
of a dnal triangular lattice: (“|) maps to dimer 

(no dimer), respectively. 

position r to a flux 17-0 at a symmetry-related position 
Xr, as shown in Fig. |l(b)[ 

X : %iv) ^ n_^iXr). (2) 

Now we consider a configuration with two fluxes inserted: 
f70(r)f7_0(Xr). Since the total flux is zero, this configu¬ 
ration can be consistently put on a finite system with a 
())-independent open boundary conditions and has a well- 
defined X-symmetry parity eigenvalue Xxi4>) = if- In 
the limit of ^ 27r, the two fluxes become two visons, 

and the ratio of parity eigenvalues Av:(27r)/Ax(0) gives 
the symmetry fractionalization of vison |29l[54||55]. If the 
symmetry X is unbroken for any (j}, the parity eigenvalue 
Xxi4>) cannot jump between -|-1 and —1, and therefore 
the vison must carry = -|-1. 

In summary, using the flux-fusion anomaly test we can 
conclude that the vison can only take a trivial symmetry 
fractionalization of = -|-1 if the flux transforms 
under X as described by Eq. ([^ if X is anti-unitary, and 
Eq. if X is unitary. 

These criteria significantly constrain the possible sym¬ 
metry fractionalizations of the vison, in the presence of 
U(l) spin rotational symmetry. However they are not 
sufficient to fix all the vison fractionalization classes. To 
be concrete, we consider the example of the kagome lat¬ 
tice, where the symmetry group is G = SO(3) xp6mxZ|’, 
with p6m being the space symmetry group of the kagome 
lattice and the group generated by time-reversal sym¬ 
metry. The fractionalization of this symmetry group on 
visons are classified by H^{G,Z 2 ) = and are labeled 
by seven Z 2 invariants: ClIi 2 , Wo-, w/, wt, and 
WctT) as listed in Table |lj Here /i, a denote two mir¬ 
ror reflections and Ip — (pa)^ denotes the plaquette- 
centered inversion, as shown in Fig. |2(a)| The first vari¬ 
able UJ 12 labels the commutation relation fractionaliza¬ 
tion [^ T 1 T 2 = ztT 2 Ti and the other six variables, in 


the form of ujx , label quantum number fractionalization 
X2 = ±1. 

TABLE I. Quantum numbers labeling different symmetry 
fractionalization of visons. The first column lists the labels 
indicating fractional quantum numbers associated with the 
algebraic relations listed in the second column. The third 
column labels symmetry operation X' used in the anomaly 
test, and the last column lists the symmetry required to de¬ 
rive X^ = -1-1. 


Label 

Algebraic relation 

X' Required symmetry 

L012 

TiTa = iTzTi 

- 

- 


= ±1 


U(l) 


cr^ = ±1 

(J 

U(l) 

LOT 

T = ±1 

T 

U(l) 

LUl 

II = ±1 

T 

Ipt. 

SO(3) 

LOfiT 

[pTf = ±1 


SO(3) 

UJaT 

{aTf = ±1 


SO(3) 


Next we divide the symmetry operations X into two 
classes according to how they act on the aforementioned 
U(l)s^ flux the symmetry operations T, and a 
satisfy either Eq. 0 or Eq. 0 (depending on whether 
the operation is unitary or anti-unitary) and can be used 
directly in the anomaly test, which immediately gives 
the constraints that vison must carry trivial symmetry 
fractionalization lot = = 1. We note that these 

symmetries satisfy t(X) = — 1 for the function t defined 
in Ref. |3S1 and the constraints X^ = -1-1 have all been 
obtained in previous works 

The other three symmetry operations Ip, fiT and aT 
do not satisfy either condition and cannot be used di¬ 
rectly in the anomaly test [they satisfy t{X) = 1]. To 
make progress, we generalize the flux-fusion argument 
by including spin rotations: we observe that if we com¬ 
bine fiT or aT with a spin rotation which reverses 

the symmetry flux (j), the new symmetry operations now 
preserve the U(l) fluxes and therefore satisfy Eq. Q. 
Similarly, satisfies Eq. ([^. Hence the combined 

symmetry X' = Xe“'®” can be used in the anomaly 
test to infer that the vison must have (Xe®^'^’^)^ = 4-1 
for X = Ip,p,T and aT. We also know that the vi¬ 
son carries = 4-1 because it has an integer 

spin, and one can further show due to the connected¬ 
ness of SO(3) group that the spin rotation must 

commute with X (see Sec. [I] of the Supplemental Ma¬ 
terial for more details). Therefore from the identity 
(Xe®'®®'^*)^ = X^(e®®®'®”)^Xe®®®'®”X“^e“®''®‘^’^ we find that 
X^ = 4-1 on visons for X = Ip, pT and aT. 

To summarize, using the anomaly test we show that in 
a SO(3) symmetric Z 2 spin liquid where the spinons carry 
a half-integer spin, the vison must carry trivial quantum 
number fractionalization except for uji 2 as listed in Ta- 
ble|Tj In fact, the symmetry fractionalization of the visons 
is completely captured by an Ising gauge theory. [The 
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Ising gauge theory is even (odd) if there are an even (odd) 
number of spin-^ per unit cell, respectively.] As listed 
in Table m deriving ojx = +1 for X = only re¬ 

quires a U(l) subgroup of the spin SO(3) symmetry, but 
deriving ujx = +1 for the other three X = Ip, aT, fiT re¬ 
quires the full SO(3) symmetry. In particular, although 
in our argument we only use the symmetry operations 
from the U(l) xi Z 2 = 0(2) subgroup generated by 
and the connetedness of S0(3) plays a crucial role 

in enforcing that commutes with X in the frac- 

tionalization class. As we will see below, the vison can 
exhibit more complicated symmetry fractionalizations if 
the spin rotational symmetry is reduced to 0(2) or a dis¬ 
crete subgroup. More details can be found in Sec. |l] of 
the Supplemental Material. 

We note that when introducing the fluxes, a branch- 
cut as shown in Fig. [^is also introduced. To consistently 
define the symmetry transformations of the fluxes, this 
branch-cut should be chosen to be symmetric under X. 
However, such a choice is impossible for a site-centered 
inversion symmetry Ig [see Fig. |2(ay| , and therefore our 
argument cannot be applied directly to fix = -|-1. This 
is discussed in more details in Sec. Hof the Supplemental 
Material. 

Z 2 SET with nontrivial vison symmetry fractionaliza- 
tion. Our argument showing the vison cannot have 
nontrivial symmetry fractionalization requires the spin 
liquid to have the full spin SO(3) symmetry. Hence 
breaking the spin rotational symmetry down to U(l) or 
0(2) = U(l) XI Z2 opens up the possibility for the visons 
to have nontrivial fractionalized quantum numbers of the 
symmetries X = Ip, fiT, aT, as listed in Table [I] Here we 
present a concrete example of Z 2 SET phases, where the 
e anyon has half U(l) charge and the vison m has both 
Ip = —1 and (/rT)^ = —1. 

We consider the Z 2 SET phase found in the BFG spin- 
1 /2 XXZ model on the kagome lattice [15]. In this model 
the spin rotational symmetry is broken from SO(3) down 
to 0(2) due to the easy-axis anisotropy. Therefore our 
argument for vison having Ip = {iTTY = (cT)^ = -|-1 no 
longer applies, although we can still use the anomaly test 
to show (/pe“^'')2 = = {aTe^^^'f = -fl. 

Indeed the vison in this Z 2 phase carries Ip = —1 and 
(/iT)2 = -I. 

To see that the vison carries Ip = —1, we notice that 
the Hamiltonian can be mapped to a quantum dimer 
model (QDM) with three dimers per site on a triangular 
lattice [46], where each spin, viewed as a hard-core bo¬ 
son, corresponds to a dimer in the latter model, as shown 
in Fig. |2(b)| In this mapping the sites of the kagome 
lattice map to the bonds of the dual triangular lattice, 
and the sites of the dual triangular lattice map to the 
center of the hexagons of the kagome lattice. There¬ 
fore the plaquette-centered inversion Ip on the kagome 
lattice [see Fig. |2(a)| becomes a site-centered inversion 
^duai triangular lattice. It is well-known that in 


the Z 2 topological phase of a quantum dimer model with 
an odd number of dimers per site, the vison excitations 
are described by an odd Ising gauge theory [56ll59] , which 
implies that they carry a nontrivial symmetry fraction¬ 
alization (1^“®'^)^ = — 1 for the site-centered inversion of 
the triangular lattice [551 El. Therefore the vison exci¬ 
tations in this Z 2 spin liquid have Ip = —1. We can also 
say that the background spinons are now located at the 
centers of the hexagons although the physical spin-l/2’s 
are at sites, which is actually evident from the mapping 
to the quantum dimer model, and the symmetry frac¬ 
tionalization Ip = —1 comes from the braiding between 
the visons and the spinon at the inversion center when 


applying Ip to the configuration in Fig. 1(b) 


In fact, our generalized flux-fusion argument is circum¬ 
vented precisely by the anticommutation between 
and Ip on visons, which is allowed for 0(2) spin symmetry 
but not for SO(3). As shown in Sec. [HI| of the Supplemen¬ 
tal Material in the fractionalization class of the vison the 
spin rotation commutes with a but anticommutes 

with /i, so it anticommutes with Ip = {^a)^. Therefore 
combined with the fact that = -1-1 on the vison, 

the anomaly test gives = —/^(e”'^”)^ which 

implies Ip = —1. 

Similarly, using the flux-fusion anomaly tes t we have 

of the 


HI 


[yiTe™^ )^ = -1-1. However as we show in Sec. 
Supplemental Material, when acting on the vison 
commutes with T but anticommutes with /r, and there¬ 
fore it anticommutes with ytT. Hence we get {yiT)'^ = —1 
for visons. 

These results have interesting physical implications. 
First of all, since in a S0(3)-symmetric spin liquid visons 
have Ip = (/iT)^ = -|-1, we have shown that the Z 2 spin 
liquid in the BFG model cannot be smoothly connected 
to the Z 2 spin liquid state found in the antiferromag¬ 
netic Heisenberg model on the kagome lattice with full 
S0(3) spin rotational symmetry, without breaking the 
space symmetry |60j . Secondly, the nontrivial fractional¬ 
ization class (ytT)'^ = — 1 for visons implies the existence 
of gapless edge states protected by both T and yi, on an 
edge that preserves the reflection /i [29] . 

Conclusion. In this work we show that in a Z 2 spin liq¬ 
uid with an unbroken S0(3) spin rotational symmetry, 
the symmetry fractionalization of visons is completely 
fixed if the spinon carries a half-integer spin. The last 
condition is automatically fulfilled in many candidate 
systems of Z 2 spin liquids studied so far [51115511551155] as 
they have an odd number of spin-^ per unit cell P[I5|. 
Our result can be directly applied to simplify the numer¬ 
ical detection [^ [^ of the symmetry fractionalization 
in Z 2 spin liquids, as only the fractional quantum num¬ 
bers of one type of spinons (either e or e) need to be 
measured, and those of m and the other spinon can then 
be determined theoretically [29115T] . 

Having determined the symmetry fractionalization of 
visons, one still has a wide variety of possible Z 2 spin 
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liquids with different symmetry fractionalization classes 
of spinons. Representative wave functions for spin liq¬ 
uids with different symmetry fractionalizations are con¬ 
structed using parton constructions including Schwinger 
boson and Abrikosov fermions [11IM31ISI], as well as 
tensor network states m- However, in all these construc¬ 
tions some spinon symmetry fractionalization classes can 
only be realized as gapless Z 2 spin liquids, and we will 
leave the question of existence of gapped wavefunctions 
for those classes to future works. 

Our argument suggests that to look for spin liquids 
with nontrivial vison symmetry fractionalization, one 
should turn to systems with reduced spin rotational sym¬ 
metry. We consider the Z 2 spin liquid in the BFG model 
as an example, where the spin rotational symmetry is 
reduced to 0(2). We find that quite interestingly the vi- 
sons carry /^ = — 1 as well as = — 1- This example 

only scratches the surface of the rich possibilities of dif¬ 
ferent SET orders in models with reduced spin rotational 
symmetry, which can be realized in systems where strong 
electron interactions interplay with strong spin-orbit cou¬ 
plings. We leave the classihcation and construction of 
symmetric spin liquids in these systems to future works. 
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Supplemental Material 


In this supplemental material we provide technical details of the discussions in the main 
text. In Sec. I we explain the mathematical details of fractionalization of spin rotational sym¬ 
metry and crystal symmetries. In Sec. II we discuss the application of flux-fusion anomaly 
test to site-centered inversions. In Sec. Ill we derive the details of symmetry fractionaliza¬ 
tion in the Balents-Fisher-Girvin (BFG) model[ 10 ]. In Sec. IV we discuss vison’s symmetry 
fractionalization in gapped 'L2 spin liquids on the square lattice. 


I. FRACTIONALIZATION OF SPIN ROTATIONAL SYMMETRY AND 
CRYSTAL SYMMETRIES: SO( 3 ) V.S. 0 ( 2 ) 

In this section we provide more details on symmetry fractionalization of the symme¬ 
try groups of spin liquids, which are a direct product of the spin rotational symmetry, 
the time-reversal symmetry and crystal symmetries. Particularly we discuss the possibility 
of fractionalizing the commutation relation between spin rotational symmetry and crys¬ 
tal symmetries, and the projective representation of spin rotational symmetry and crystal 
symmetries if the former is reduced from the full SO( 3 ) symmetry. 

First, we consider the simple case where the spin liquid has a full S 0 ( 3 ) spin rotational 
symmetry. In this case it is well known that the projective representation of the spin rota¬ 
tional symmetry is classihed by Ff^(S 0 ( 3 ), Z12) = Z2, which consists the trivial integer spin 
representations and the nontrivial half-integer representations. If we enlarge the symmetry 
group to S 0 ( 3 ) X Z2, where Z2 = { 1 ,X} represents the time-reversal symmetry or a crystal 
symmetry, its fractionalization classes are given by the Kunneth formula 

H^{G X Z2,Z2) = H\G,Z2) X Z 2 X H\G,H\Z2,Z2)) 

( 1 ) 

= H\G,Z2) XZ2 X ipi(G',Z2), 

where G = S 0 ( 3 ) and we have used Ff^(Z2,Z2) = Z2 and if^(Z2,Z2) = Z2. In other words, 
the fractionalization of the symmetry group SO( 3 ) x Z2 are labeled by the combination of 
the fractionalization of G, the fractionalization of of Z2: if^(Z2,Z2) = Z2, and an extra 
factor if^(G,Z2), which classihes group homomorphisms from G to Z2. For a continuous, 
connected group G, the image of any g & G under the homomorphism can be continuously 
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connected back to the image of the identity. Since now the images are valued in a discrete 
group Z 2 , we immediately conclude that the homomorphism can only be trivial (since the 
image of the identity must be identity). More intuitively, we observe that last factor can 
be thought as the fractionalization of commutation relation between an element g of G and 
the nontrivial element X of Z 2 , f{g) = gXg~^X~^, and f{g) takes value in ±1. Since the 
values of f{g) are quantized, we must have f{gi) = /(fi' 2 ) if two elements gi and g 2 can 
be smoothly connected in G. Particularly if G G is connected to identity we conclude 
fid) = /(I) = +1- Since G = SO(3) is connected, all its elements must commute with X 
in any fractionalization class. We notice that the same argument was given in Ref. [9]. The 
conclusion is used in the main text to show that = +1 and = +1 implies 

= + 1 . 

The symmetry fractionalization become more complex if the spin rotational symmetry is 
reduced to U(l) x Z 2 = 0(2). First, the fractionalization of the spin rotation itself become 
more complex as iP^( 0 ( 2 ),Z 2 ) = Z^, and are labeled by two independent Z 2 variables 
denoting = ±1 and (e*^'®^^)^ = ±1, respectively. So if we still assume that the 

e particle carries a half charge of the U(l) symmetry = —1, the vison can take 

(gi7rSJ')2 _ regardless of whether the e particle takes (e*’^'^^)^ = ±1. But with time- 
reversal symmetry, we can exclude the possibility that both e and m have (e*^"^*^)^ = — 1. 
Second, the spin rotation can anticommute with a crystal symmetry X because it is 
no longer connected to the identity in 0(2). Therefore if the spin rotational symmetry is 
reduced to 0(2), our argument showing = -|-1 from = -|-1 no longer holds and 

the vison can take nontrivial symmetry fractionalization. 


II. SITE-CENTERED INVERSION. 

In this section we explain that our argument of vison carrying trivial inversion symmetry 
fractionalization uj = +1 only applies to a plaquette-centered inversion Jp, not to a site- 
centered inversion Ig, for lattice models with a half-integer spin on each site [1]. The reason 
that our argument fails for site-centered inversion G is that the branch-cut connecting the 
two fluxes shown in Fig. 1(b) of the main text cannot be chosen to be symmetric under 
Is- By dehnition the branch-cut separates sites where the unitary transformation creating 
the fluxes acts discontinuously and thus cannot go through a site, however without going 
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though the site at the inversion center it cannot be inversion-symmetric. 

To £x this issue, we can choose a branch-cut that is not invariant under inversion, and 
after the inversion symmetry operation we apply an unitary transformation to restore the 
original position of the branch-cut. Particularly if we choose a branch-cut such that it and 
its inversion image encloses only the site of the inversion center, we can apply the unitary 
transformation at the inversion center to restore the branch-cut after the inversion. 

Therefore the anomaly test works for the combination of inversion symmetry Ig (plus 
and this unitary transformation. If the site at the inversion center carries an integer spin, this 
unitary transformation becomes trivial at 0 = 27r and we still get the result that If = -|-1. 
However, if the site at the inversion center carries a half-integer spin, at 0 = 27r this unitary 
transformation gives an extra —1 Berry phase. This explains that in such models vison 
always carries opposite fractional quantum numbers [2], and our argument fixes 

them at = -|-1 and I] = —1, respectively. The relation I] = —can also be derived 
directly from odd Ising gauge theories and vison’s PSG on different lattices [3-8]. 

Since P represents a 27r rotation, = — 1 for the vison can be interpreted as the vison 
braiding around a background spinon (e or 0) sitting at the inversion center, and = —If 
is fixed by the spin-1/2 per unit cell. The anomaly test shows that for SO(3) symmetry, 
background spinons must sit at the sites. When the spin symmetry is reduced, the position 
of the background spinons can shift away from the sites to the plaquette centers, as we see 
in the example of the BFG model discussed in the main text. 


III. SYMMETRY FRACTIONALIZATION IN THE BFG MODEL 

In this section we systematically calculate the symmetry fractionalization of the vison in 
the BFG model [10]. The results obtained here are used in the main text to argue that in 
this Zi 2 spin liquid state the vison carries nontrivial symmetry fractionalization class with 
Ip = —1 and {prY = —1- The symmetry group of the BFG model is 0(2) x xp6m, where 
the spin rotational symmetry 0(2) = U(l) x Z 2 is generated by Sz and respectively. 

As shown in Table I of the main text, the fractionalization of translation symmetries are 
labeled by U 12 and it is separated from all other symmetries, so we only need to consider 
the point group symmetry generated by fi and a. 

We start with the spin rotational symmetry and the point group symmetry. We note 
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that as an on-site nnitary Z 2 symmetry, it is not easy to obtain the fractional qnantnm 
nnmber of = ±1 directly from gronnd state wave fnnctions [11, 12]. Hence here 

we take an indirect approach and calculate this fractional quantum number using algebraic 
relations for the fractionalization of and the other point group operations for visons. 
For simplicity, we use ux to denote fractional quantum numbers: oogy = 
uj„ = ui = Ip, Ugy^ = i^sya = and Ugyi = respectively; 

we use axY to denote the commutation relation fractionalization between two symmetries 
X and Y: a = and Ip = agyilpC^'^^^ . Using the 

algebraic relation (XY)'^ = we have 

^syX ^sy^syX 1 ^ h) ^P" (^) 

Furthermore, using the relation Ip = (a/i)^ we derive the following relation, 

^ syl ^ syy^ sya • ( 3 ) 

Next, we £x the fractional quantum number of crystal symmetries cn^, uj, ojgyp, 
Ugya and Ugyi using the symmetry eigenvalues of two-vison wave functions. In this exactly 
solvable model visons are located in the triangles of the kagome lattice, and a wave function 
containing two visons located in two triangles i and j is obtained by applying the following 
string operator on the ground state [10], 

n 25", (4) 

r£p{ij) 

where the product is taken along a path p{ij) connecting triangles i to j, as shown in 
Fig. 1. The path can only go straight or make “±30°” turns, and the ± sign in the front 
represents a Z 2 gauge choice depending on the path p\ if two paths enclose N hexagons, the 
corresponding string operators differ by a sign (—1)^ because of the constraint that there 
are three up spins and three down spins on each hexagon [see Fig. 1(c)]. In this model 
there are two types of visons, denoted as red and blue visons in Ref. 10, located in the up 
triangles and down triangles of the kagome lattice, respectively. It is easy to check that a 
string connecting two visons with the same(opposite) colors always contains an even(odd) 
number of S^. 

We begin with the mirror reflection a and consider a string operator along the path 
p{ij) connecting two triangles j = o'(f), as shown in Fig. 1(a). The string operator itself is 
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a 





FIG. 1. String operator creating two visons. The dashed line shows a path p{ij) connecting two 
visons, marked by cross symbols in two triangles of the kagome lattice, where the red and blue 
color labels visons in up and down triangles, respectively. The string operator is a product of 
operators on the sites marked by black dots, as in Eq. (4). (a) A vison string operator symmetric 
under mirror reflection a. (b) A vison string symmetric under mirror reflection p. (c) A string 
operator connecting two visons symmetric under the 180-degree rotation Ip. After Ip the string 
in black color is mapped to another string in grey color, and the two strings differ by a factor — 1 
because together they enclose one hexagon. 


symmetric under a, implying that Uo- = +1. Because a maps up triangles to up triangles, 
the string operator connects two visons with the same color and contains even number of 
S^, and it commutes with [see Fig. 1(a)]. Therefore the two-vison wave function is 

symmetric under implying that cOgya = +1- 

Similarly for the mirror reflection /r, we consider a string operator along p{ij) where 
j = as shown in Fig. 1(b). The string operator is still symmetric under /r, hence 

Up = +1. However since p maps an up triangle to a down triangle and vice versa, the string 
operator connects two visons with opposite colors and contains an odd number of and 
therefore anticommutes with This implies that Ugyp = — 1. 

Lastly we consider the inversion symmetry Ip. In this case the path p{ij) is mapped to 
another path I[p{ij)], and together with p{ij) it forms an inversion-symmetric loop that 
encloses an odd number of hexagons because Ip is centered at a hexagon, as shown in 
Fig. 1(c). Hence the two paths differ by a minus sign, implying that the two-vison wave 
function is antisymmetric under Ip and therefore ui = —1. Furthermore the length of this 
string is also odd because it connects a up triangle and a down triangle, so the wave function 
is also antisymmetric under . Combining these two results we see that the two-vison 
wave function is symmetric under Ip and therefore Ugyi = +1. 
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So far using the vison string operators we explicitly derive the following symmetry frac- 
tionalization = +1, Ugy^ = —1, cosya = +1, = —1 and Usyi = +1. We note 

that the results of uji and Ugyi are consistent with the results we obtained using 

general arguments in the main text. Now substituting the results of uj„, Ugyy and Ugya 
into Eq. (2), we get 

^sy^syy f? ^sy^sya “1“1. (h) 

Multiplying these two equations together, we see that Ugy is canceled and using Eq. (3) we 
get asyydsya = o'syi = —1. Finally putting this result and ui = —1, Ugyi = +1 into Eq. (2) 
we get Usy = +1. 

These results can be understood intuitively by associating the color of vison with its sym¬ 
metry quantum number = ±1. Since Ugy = +1, the visons carry linear representations 
of and the vison on each triangle carries one of the two possibilities = ±1. Then 
we consider the string operator in Eq. (4) that creates two visons. Because an even-length 
string creates two visons with the same color, and the even-length string is invariant under 
the global symmetry transformation two visons with the same color have the same 

quantum number of . Similarly as an odd-length string creates two visons with oppo¬ 
site colors and it changes sign under the red visons and the blue visons have opposite 

giTrSJ' quantum numbers. Without losing generality we assume a red (blue) vison carries 
gjTrSJ' _ _j_]^ (—1), respectively, (the two choices are related by a Z 2 gauge transformation 
and thus physically indistinguishable). This assignment of quantum number is consistent 
with the commutation relation fractionalization (Jsyi = — 1, since Ip maps a red vison to a 
blue vison and therefore anticommutes with This also implies that asyp = — 1 and 

o'sya = +1, because fx flips vison’s color and a preserves vison’s color. These are consistent 
with what we can get from Eq. (5) and cogy = +1- 

In summary, we have shown that the visons in the Z 2 topological state of the BEG 
model [10] have symmetry fractionalization = -|-1 and anticommuting with Ip. 

This explains how the result of {IpC^'^^^Y = +1 from the anomaly test is consistent with 

P = - 1 . 

P 

Next we consider symmetry fractionalization involving the time-reversal symmetry. We 
first notice that when the BEG model is mapped to a quantum dimer model (QDM), the 
time-reversal symmetry T in the spin model is mapped to a rather unconventional time- 
reversal symmetry in the dimer model(i.e. it flips the dimer occupations). In fact, T = 
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T^^T^sy ^ which is just the complex conjugation K, is the more natural time-reversal operation 
in the dimer model which do not change the dimer conhgurations. Using this correspondence 
the symmetry fractionalization in BFG model can be derived from that of the QDM. 

First, it is well-known that in the QDM the vison has = -|-1, thus in the BFG model 

we get = -|-1. Furthermore we have = -|-1 from the flux-fusion anomaly test, 

and = -|-1 from the previous discussions, so we conclude that T commutes with 

^insy 

Second, in the QDM the vison takes trivial symmetry fractionalization {f^TY = +1 and 
{aTY = +1, which map back to = +1 {aTe'^^^^Y = +1 ^he BFG model, 

consistent with the results of the generalized flux-fusion anomaly test used in the main text. 
On the other hand, since commutes with T and a but anticommutes with p, we get 
that commutes with aT but anticommutes with fxT. So we hnd (uT)^ = -|-1 but 

(/iT)2 = -1. 

IV. SYMMETRY FRACTIONALIZATION ON A SQUARE LATTICE 

In this section we discuss the symmetry fractionalization for visons on a square lattice, 
whose space group is pAmm. We will also assume there is a half-integer spin per unit cell, 
and call the bosonic spinon e. The algebraic relations and the corresponding invariants are 
dehned in Table I (here we largely follow the notations in Ref. 9 and 13, but we use bond- 
centered mirror reflection and plaquette-centered inversion I = PxPxy to ensure that in 
the flux-fusion anomaly test a symmetric branch-cut can be drawn without going through 
a site, as discussed in Sec. II). 

With U(l) 5 ^ symmetry, we can immediately £x 

~ ^Pxy ~ = 1 ; = ~ 1 - ( 6 ) 

It is shown in Ref. [14] that = 1. 

We now consider the other five invariants, u/, Ctyp^-, o'Tpx,-, ^Tpxyi ^Ttx- To fix those we need 
to consider SO(3) spin symmetry. Using the same argument as in the main text, we fix 

= (Ji = 1. 

To fix arty (which is equal to we put the system on a cylinder geometry where the 

y direction is compactihed, and assume the circumference Ly = 2 (mod 4). We will consider 



TABLE I. Algebraic relations on square lattice. The first seven can be constrained by the same 
argument for the case of triangular lattice. The last three are only present in square lattice. 


Algebraic relation 

Invariants for visons 

Required symmetry 

T T - hi 

(^txty 

U(l) 

= ±1 


U(l) 

= ±1 

^Pxy 

U(l) 

{PxPxyY = Tf 

cri 

SO(3) 

T—1 

-H 

II 

U 

CJj' 

U(l) 

TP^T-^P^ = ±1 

^Tp, 

SO(3) 

TPxyT~^Pxy = ±1 

^Tpxy 

SO(3) 

TxPxTxPp^ = ±1 


U(l) 

TxPyTx ^Py ^ — T1 

'^txPy 

U(l) XI ^2 

TTxT-^T-^ = ±1 


U(l) 


the invariant 

_ rj^2 _ ^rjnLy/2Y _ rjnrjnLy/2rjn-l ^rjnLy/2^-l ^ 

All the quantities are computed on Schmidt states for an entanglement cut along y. We view 
this result as an invariant of the dimensionally reduced ID SPT, where now Ty is viewed as 
an on-site symmetry in the ID system. Using the flux-fusion anomaly test, we can easily 
prove that = 1 for the vison sector. We have already known that = 1 for the 

vison. is related to the momentum polarization, which is then determined by the 

topological twist factor of the anyon sector. In general is not a ID SPT invariant, 

but as shown in Ref. [12], in the presence of Py it is quantized exactly to the topological 
twist factor (as long as the magnitude does not vanish), which is 1 for the vison. So we 
conclude that = 1 for the vison sector. Since we have chosen 

a cylinder with Ty/2 being an odd integer, we have shown that (JTty = 1- 

A similar argument works for Ut^py. We still assume Ly is even and consider the translation 
along X. Dehne c^t^py as the ratio of the eigenvalues of Py = PyC™^^ per unit length in the 
\m) and |1) sectors. Aagain U(l) 5 ^ flux threading hxes Ut^^py = 1, and since we assume Ly 
is even, the eigenvalue of per unit length is 1. We therefore conclude that at^p^ = 1. 
Notice that in this argument we do not need to use the connectedness of SO(3), so it would 
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work equally well if we break the spin rotational symmetry down to U(l) xi Z 12 . 

In summary, with SO(3) spin rotational symmetry and a half-integer spin per site the 
symmetry fractionalizations of visons are completely hxed: except = —1, all others 
must be trivial. 
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